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Abstract 

We study in an effective operator approach how the effects of new physics from various scenarios 
that contain an extra Z' neutral gauge boson or doubly charged scalars, can affect and thus 
be tested by the precision polarized Moller scattering experiments. We give Wilson coefficients 
for various classes of generic models, and we deduce constraints on the parameter space of the 
relevant coupling constants or mixing angles from the results of the SLAC E158 experiment where 
applicable. We give also constraints projected from the upcoming 1 ppb JLAB experiment. In the 
scenario where the extra Z' is light {Mz 1 <C M^), we obtain further constraints on the parameter 
space using the BNL g — 2 result where it is useful. We find that the BNL deviation from the 
Standard Model cannot be attributed to a light extra Z' neutral gauge boson. 
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I. INTRODUCTION 



Polarized electron scattering has a long and illustrious history that began with the cel- 
ebrated SLAC polarized eD scattering [1], which was instrumental in putting the Sandard 
Model (SM) on a firm experimental ground. The latest experiment of this class is the polar- 
ized M0ller scattering. As only the electrons are involved, theory calculations at quantum 
loop levels can be performed free of hadronic uncertainties making it an ideal process for 
precision SM testing; any deviation from expectations will also be a clear signal of new 
physics beyond the SM. 

A while ago, an experiment of this type has been successfully completed [2], where parity- 
violating measurements have been made from the scattering of longitudinally polarized elec- 
trons - either left-handed (LH) or right-handed (RH) - on an unpolarized electron target. 
At a momentum transfer of Q 2 = 0.026 GeV 2 , the left-right asymmetry 

_ da L - da R 
daL + da R 

has been measured to an accuracy of approximately 20 ppb, and the weak angle determined 
to be sin 2 6 W = 0.2397 ± 0.0010 (stat.) ± 0.0008 (syst.), with its running established at the 
6a level. In the SM, A^ R arises from the interference between the electromagnetic and the 
weak neutral currents, and is given by 



4 sm _ G »Q 2 1 - y 



l-4sin 2 ^ s 



Q 2 = ys, s = 2m e E heam , (2) 



LR V2nal + y 4 + (l-yy 
where = 1.16637(1) x 10 -5 GeV -2 , a -1 = 137.0356, and terms involving the electron 



mass m e have been dropped when appropriate. T 



experimental accuracy and more is given in Ref 



.a. 



ie electroweak corrections matching this 
The value of the weak angle at Q 2 = 
4| . We use this value in 



is found to be sin 2 9$ s = 0.23867 ± 0.00016 in an updated analysis 
all our numerical calculations below. 

Recently, a 12 GeV energy upgrade to the intense electron facility at JLAB has been 
approved. This raises the possibility of measuring Alr to an accuracy of 1 ppb. Since the 
theoretical uncertainty can in principle be calculated to this accuracy, a successful measure- 
ment will be sensitive to the new physics at the several TeV scale, which is also the physics 
landscape to be explored by the LHC. 

Also recently, an anomalous abundance of positron flux in cosmic radiation has been 

n 

reported by the PAMELA collaboration [5J. This confirms the excess previously observed 



by the HEAT and AMS Q experiments. The ATIC 8] and PPB-BETS |e| balloon 
experiments have also recorded similar excess. One possible origin is from a hidden sector 
and low scale dark matter annihilation 10j. The hidden sector scale is O(GeV), and so is 
at the opposite end of the new heavy physics. The characteristic of this class of model is 
that possible couplings to SM matter fields must be extremely weak. A very simple model 
was constructed previously to study this possibility [111 ] . and M0ller scattering was found to 
give very stringent constraints. 

Motivated by all these encouraging developments, we examine in detail how an im- 
proved Alr measurement could shed light on new physics. We focus here on classes of 
non-superymmetric new physics and how they might manifest themselves in precision po- 
larized M0ller scattering. It is well known that A^r can be used to probe the existence 
of doubly charged scalars and extra Z' bosons. The doubly charged scalars are interesting 
in their own right, and they also arise in Type II seesaw models for neutrino masses 12]. 
For extra Z' bosons, the literature is particularly rich with well motivated models. These 
include models of grand unified theory (GUT) based on the E 6 or SO(10) gauge group |l3| . 
hidden sector Z' models fill. [l4 ] . and left-right symmetric [l^ and warped extra-dimensional 
Randall- Sundrum (RS) |l6[] models based on the gauge group SU(2) L x SU(2) R x U(1)b-l- 
The effects of a new particle such as the doubly charged scalar of the extra Z' boson can 
be encoded in an effective operator approach, which for Am involve only two dimension six 
operators: 

1 = erfleer (l Xe, 2 = erfRee^Re , L = -(l- 75 ), fl = I(l + 75 ). (3) 

In the next section, we describe this effective operator approach. 

The rest of the paper is organized as follows. In Sec. 1111} we discuss the classes of new 
physics models. Besides updating relevant known results, we also present a new calculation 
of the contributions in the Minimal Custodial RS (MCRS) model with bulk symmetry 
SU(2) L x SU(2) R x U(1)b-l, and we give comparisons with the usual four dimensional 
left-right symmetric models. Sec. [IV] contains our conclusions. The details of the MCRS 
model and the specifics of the fermion localization are collected into Appendix. [A] 
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II. EFFECTIVE OPERATORS 



A. High Scale New Physics 

If the new physics is at a scale higher than the Fermi scale v ~ 247 GeV, an effective 
operator approach is powerful and avoids much model dependence. Denoting the scale of 
new physics by A and assuming that the SM gauge group holds between A and v, the effective 
Lagrangian with dimension six four-lepton operators added is given by 

£ = C-SM + ^6 , 

+ ^(L a ^L b )(L b ^L a ) + h.c, (4) 

where a, b are SU (2) indices, and L = (z/, e) \. The subscripts L and R denote the usual 
chiral projection. We have included only terms relevant for the calculation of A^ R , i.e. the 
first generation leptons. Note that for scattering of a polarized beam on a polarized target, 
there is also a scalar operator that contributes: 

C s = ^Le R e R L + h.c. (5) 



A full analysis of dimension six leptonic operators can be found in Ref. [17] . 

The Wilson coefficients c and d are given by specific models at A. One has to do a 
renormalization group (RG) running analysis to determine their values at the relevant low 
energy scale where the experiments are performed. The RG equations for these coefficients 
can be found in Ref. Numerically, we find that the RG running effects are not significant. 
They give about a 4% effect running between A (typically a few TeV) and v , and hence can 
be neglected in the first approximation. Below v the running is governed by the well known 
Pqed funtions. Including them leads to a change of about 10% in the values of the coefficients 
between A and a/s ~ 0.1 GeV. We neglect these small RG running effects below which would 
not impact our analysis given the current level of precision in the experimental inputs and 
the exploratory nature of our study. 

The asymmetry A LR is found by calculating the diagrams depicted in Fig. [1] The shaded 
blob schematically includes all new physics effects such as Z-Z' mixing and virtual exchanges 



4 



FIG. 1: Feynman diagrams contributing to Alr- The shaded blob contains all new physics effects. 



of new particles. For A ^> s ^> m e valid for the 12 GeV experiment, one has 



A LR = A™ + 8 A LR 

4G> y(l-y) 



1 c - 2 fl MsA i C LL C RR 

— — sin y w H = 

4 w J 4V2G U A 2 



(6) 



where y = — c^ L = c^l + g?ll, and 5^^^ denotes the deviation of the asymmetry from the 
SM prediction arising from new physics. The one-loop SM radiative corrections are encoded 
in the precise value of sin 2 #\y evaluated in the MS scheme. The effective operators, (9i,2, 
can arise from numerous new physics models, and those coming from the tree-level exchange 
of virtual particles are the most important. We thus examine models with at least an extra 
neutral gauge boson, Z' ', or doubly charged scalars below in Sec. IHIl where we give Wilson 
coefficients in various models. Of course there could be extra charged gauged bosons that 
also contribute in specific models, but these typically come in at the loop level and so are 
less important. 

We remark here that the new physics effects in M0ller scattering can also come in through 
the Zee coupling being modified. This is particularly evident in cases when the heavy states 
of the new physics also have very small couplings to electrons so that the virtual exchange 
process is suppressed both by the propagator effect and the small coupling. Example of 
this can be found in the extra hidden sector/shadow U(l) and RS models which we discuss 
below in Sec. IfflDl and lITTEl 



B. Low Scale New Phyiscs 

For the class of models where new physics arise at a scale A < v, two effects are at work. 
The first is due to the exchange of the new low mass particle. The propagator effect has to 



be compensated by its small couplings to the electrons of different chiralities. The second 
is due to the modification of the SM Z coupling to the electrons. Well known examples are 
models with extra Z' that can mix with the Z. Details are given below in Sec. IHIl 



III. MODEL CONSIDERATIONS 

In this section, we discuss how new physics can be probed by polarized M0ller scattering 
in various non-supersymmetric extensions of the SM. We begin with models that contain 
extra neutral Z' bosons. A comprehensive review on the physics of Z' can be found in 
Ref. 



Irrespective of the ultimate origin of the extra neutral gauge boson, its interactions can 
be described by the gauge group, SU{2) L x U{1)y x U(l)', at a scale not too high compare 
to v [12(1 . Denoting the extra U(l)' gauge field by X^, its interaction with the SM matter 
content is given by 

£ = -\x„X^ - e -B, v X^ + g'zff-ffX^ - \m\X^ + \D,H\ 2 , (7) 

with the covariant derivative given by 

D, = d,- ig L T a W; - i^B, - i 9 -z H X, , (8) 

where gi, gy, and g' are the gauge couplings of the SU(2) L , U(1)y, and U(l)' respectively, 
/ a species label that runs over all SM fermions, and Zf (zh) the SM fermion (Higgs) charge 
under the U(l)'. Other notations are standard. The paramter e characterizes the kinetic 
mixing between B^, the hypercharge gauge field, and X^, and is expected in general to be 
of order 10" 4 to 10" 2 fl Q. 

The kinetic energy terms for the gauge bosons can be recast into canonical form through 
a GL(2) transformation 

' (9) 





e 



2 



In most cases involving high scale new physics, the kinetic mixing is negligible. However, 
this is not true in general for low scale cases. The mass term in Eq. (J7D can arise from a 
hidden sector scalar interactions that breaks the U(l)' gauge symmetry. The details are not 
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important for this study. After the spontaneous breaking of the electroweak symmetry, this 
mass term will lead to a mass mixing among W3, B, and X fields which depend on zh- 

For models with zh = but a non-vanishing kinetic mixing between the neutral gauge 
bosons, the transformation between the weak and mass basis is given by 



w 3 
\ x '/ 



cw —sw 0^ 
sw cw 
I) 



1 

Cri 

s n 



\ 



"-"n J 



/ 7 \ 

Z 

\ z '/ 



(10) 



where sw {cw) denotes sm9 w (cos9 w ), and similarly for the angle rj. The first rotation is 
the standard one that give rise to the would be SM Z, which we label as Z , while the second 
diagonalizes the mixing between Zq and X' that produces the two neutral gauge eigenstates, 
Z and Z' . This mixing angle is given by 

2s w s t 



tan 2r] 



x 2 + s w s 2 e 



x 



cwM. 
Mi 



x 



w 



(11) 



and the masses of Z and Z' are given by 



Ml 



x 2 + 1 + s 2 w s 2 =f 



x 2 - 1 + s 2 w s 2 ) 2 + As 2 w s 2 



(12) 



where M_ (M + ) denotes the lighter (heavier) of the two states. For high scale new physics, 
My/ <C Mx- Since s e < 10~ 2 , the mixing angle 77 is very small as rj < e when s e is small and 
x > 1. For M^y ^> Mx, e and 77 can be of the same order. The relative sign between them 
depends on the relative magnitude between x 2 and s 2 . 

Models in which e is vanishing but zh is not are more typical of high scale new physics 
scenarios. Since there is now no kinetic mixing, the Zq mixes with the X. The mass matrix 
of this system can be written as 



M 2 Zo _ x 



Mf o A 2 
A 2 M 2 



(13) 



with the mixing angle and the masses of the resulting Z and Z' given by [18|, [l 

2A 2 



Mf -M 2 ' 



and 



M 



z, z 1 



tan 29 



M 2 Zo + M 2 xT \l{M 2 Zo -M 2 x ) 2 + ^ 



(14) 



(15) 



Note that both A and Mx are model dependent. For example, if the SM Higgs is charged 



under both U(l) factors, A 2 = — \zhV % g'\J 9 2 + 9y an d M x = \v 2 g' 2 z 2 H + Mq where M is 
the mass from higher scale. 

In either case, the Wilson coefficients can now be calculate from e, g', Zf, zh, and A in 
a given model. For quick reference, we summarize here the results from below the Wilson 
coefficients for high scale new physics in Table [H 

New Physics c' LL c RR 

3ttu I „ i / 5 _ \ 3ira / 1„, / 5 



e 6 gut fg(^ + V^J ^[*<*-y/v°f> 

rr/,\ v TT (-\\ 1 9y 1 (29y-9r) 2 

U(1)r X U(1) B - L 4 

U(1) B -l g' 2 g' 2 

U(l)x g' 2 z 2 L g' 2 z 2 

P±± -\Vee\ 2 /2 



TABLE I: Summary of Wilson coefficients c' LL and crr for various generic high scale new physics 
models containing an extra Z' neutral gauge boson or doubly charged scalars (P ±± ). 



A. GUT Models 

The appearance of extra U(l) gauge groups is exemplified by the E$ GUT model. One 
way they can arise is through the symmetry breaking chain 

E 6 -> 50(10) x U(l)ip -> SU(5) x U(l) x x U{1)^ . (16) 

The £7(1) factors here are not anomalous. There are also no kinetic mixing, and the mass 
mixings between the gauge bosons of the extra U(l) factors and the SM Z are phenomeno- 
logically found to be negligible. The X boson above is now a linear combination of the ip 
and x gauge bosons 

X»(P)=B»fanP + B£cosP, (17) 
where B^ are gauge fields of the U(1)^ )X respectively. The X M is taken to be _ the lighter 



of the two neutral linear combinations. From the fermion charges given in Ref . 18|, [2JJ , we 
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find 



c LL - c RR 



Ana 

—c p | Cp + 



3c 



w 



(18) 



where cp (sp) denotes cos/? (sin/3). 

Denning 5 LR = (A LR - Af%)/Af% = 5A LR /A S L ^, we plot in Fig. [2] the resulting 5 LR due 
to GUT new physics. We plot also the constraints from the Tevatron new physics searches, 



E158 














0.5 1.0 


\\) ^ ^y/ -0.5 





FIG. 2: The ratio, 5lr, as a function of sp. The labelled solid lines are the SLAC E158 upper 
limit, Slr = 0.122, and cases for CDF lower mass limits on Zj and Z', A = 729, 891 GeV. The 
other solid, dashed, and dotted lines denote cases where A = 1, 2, 4 TeV respectively. The narrow 
band marks the projected JLAB constriant, \5lr\ < 0.0136. All limits are given at 95% CL. 



the SLAC E158 experiment, and that projected from the the JLAB 1 ppb measurement. We 
use the Tevatron limits on the Z\ and Z' v masses to illustrate the current collider constraints 
on the new physics scale A for the GUT extra neutral gauge bosons 2 . The 95% confidence 
level (CL) lower mass limit for Z\ (ZQ is M z / > 729 GeV {M z < n > 891 GeV), representing 
the lowest (highest) exclusion limit in the CDF search 22]. The SLAC E158 measurement 

n 

of the left-right asymmetry [2| gives 



A LR = 131 ± 14 (stat) ± 10 (syst) ppb . 



(19) 



1 As can be seen from Eq. $HD, with A - 0(100) GeV and M x > 0(1) TeV expected, the Z - X mass 
mixing is small with sin6> < 0(0.01). We thus follow common practice and neglect the small corrections 
to Alh due to the mass mixing and the concurrent modification in the gauge-fermion couplings, which 
are proportional to (sin6*) 2 . 

2 In our notation, Zj and Z' v are just X with tan/3 = y/3/5 and — -y/5/3 respectively. See Eq. (fT"7|) . 
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Taking Aff = 1.47 x 10~ 7 for Q 2 = 0.026 GeV 2 and y = 0.6, one can infer from this that 
—0.337 < 5lr < 0.122 at 95% CL 3 . Lastly, as a benchmark assume that the total la error 
of the new JLAB measurement of A^ R to be 1 ppb with a central value given by the SM 
value, Afj^, one expects in such a case \5l R \ < 0.0136 at 95% CL. We see from Fig. [2] that 
current experiments do not constrain TeV scale GUT new physics, but the upcoming 1 ppb 
JLAB experiment can be expected to do so. 



B. Left-Right Symmetric Models 

For left-right symmetric models, the gauge group is taken to be SU{2) L x SU{2) R x 
U(1)b-l- The SU(2) R is broken down to a U(1) R by a triplet Higgs field at a scale greater 
than a few TeV, leaving SU(2) L x U(1) R x U(1) B -l- In this case, the U(1) R x U(1) b ~l is 
equivalent to the U{l)y x U(l)' [3]. The charges Zf and z B can be determined by enforcing 
anomaly cancellation conditions (see below). We find 

c' L L-CR R = \{3g 2 Y -g R ) , (20) 

where g R is the gauge coupling of SU(2) R . Interestingly, note that even if g R = g L , 5l R will 
not vanish. We plot in Fig. [3] the resulting 5lr for left-right symmetric models, as well as 
the relevant SLAC E158 constraint and that projected from the JLAB 1 ppb measurement. 
As a comparion, we plot also the curve that arises from the Z^ R 95% CL lower mass limit, 
Mz LR > 630 GeV, obtained from the combination of electroweak data, direct Tevatron and 



indirect LEP II searches 



1< 



| . We see that current experiments do offer, if only barely, some 
constraints on the TeV scale new physics here compared to the GUT case. But for useful 
ones, the precision of the upcoming JLAB experiment is still much needed. 



C. SU(2) L x U(1) Y x U(l) x Models 

There are interesting bottom up models where the extra U(l) factors are not immediately 
related to the remnant of some broken down non-abelian groups. The prime example is the 
SU{2) L x U{1)y x U{1)b-l model. The assumption is that the SM is charged under U(l)x, 

3 We have followed the more common practice here in taking the ler total error as the statistical and 
systematic uncertainties added in quadrature. The 95% CL corresponds to an error interval of ±1.96er. 
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FIG. 3: The ratio, Slr, as a function of g^. The labelled solid lines are the SLAC E158 upper 
limit, 5lr = 0.122, and the case for the Zlr lower mass limit, A = 630 GeV. The other solid, 
dashed, and dotted lines denote cases where A = 1, 2, 4 TeV respectively. The narrow band marks 
the projected JLAB constraint, < 0.0136. All limits are given at 95% CL. 

and since this group is known to be anomalous, one has to demand anomaly cancellations for 
chiral fermions under the full SU(3)c x SU(2)i x U(l)y x U(l)x- As neutrinos are massive, 
we introduce one sterile neutrino, Nr, per family and denote its charge by Zn- Assuming 
these charges are family independent, we have seven independent charges Zf, z x , and z H . 
The anomaly conditions arise from 

[SU(3)] 2 U(l)x, [SU(2)] 2 U(1) X , [U(1) Y ] 2 U(1) X , [U(l) x ] 3 . (21) 

The mix gravitational and U(l)x anomaly give the same condition as [?7(l)x] 3 and is re- 
dundant. The Yukawa terms that give quark and lepton masses yielded one additional 
condition. Thus there are two independent charges, which we choose to be zl and z e . The 
other charges are 

12 4 
zq = ~^ Zl ' Zu = 3 Zl ~ Ze ' Zd = ~~3~ Zl ~^ Ze ' Zn = ^ Zl ~ Ze ' z h = z l — z e ■ (22) 

For the interesting case where X = B — L, we have the familiar lepton and baryon 
number assignments, and Zh = 0. More importantly, if the extra Z' arises directly from the 
breaking of a gauged U(1)b-l, M0ller scattering will yield only the SM result, since in that 
case zl = z e \ this leads immediately to c' LL — crr = 0. In general, zl ^ z e , and we find 

c' LL - crr = g'\z\ - zl) , (23) 
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where g' is the gauge coupling constant of the U(l)x- Note that in this analysis the masses 
of the sterile neutrinos do not come into play. Hence the details of the mechanism for 
generating active neutrinos masses will not change our results. We plot in Fig. H] contours 
of constant 5lr as a function of g'zi and g'z e for A = 1 TeV. For other values of A, the 
contour plot is simply rescaled by a factor of (A/1 TeV) as can be seen from Eq. 
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FIG. 4: Contours of constant 5lr for A = 1 TeV. Note the contours for the SLAC E158 lower 
and upper limits, 5lr = —0.337 and 5lr = 0.122, and that for the projected JLAB limits, 5lr = 
±0.0136. All limits are given at 95% CL. 



D. Hidden Sector/Shadow Z' Models 

Another class of models with an extra U(l) factor, which has received renewed attention 
recently due to the observation of excess cosmic positron flux, is one such that there are no 
direct couplings of the extra gauge field to SM fermions or Higgs fields. The extra £7(1)/^ 
is a gauge symmetry of a hidden sector with its own scalar sector. The hidden sector and 
the SM communicate through hidden scalars couplings with the SM Higgs [llj and kinetic 



mixing 



23]. In this class of models zh = 0, and the fermion couplings to Z and Z' are given 
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by [ni 




(25) 



(24) 



where Q L z ,R {f) = T z (f LyR ) - s 2 w Q f with Q f = Tf + Y f /2. We see that the neutral current 
couplings are not only rotated as indicated by the c v factor, but they also contain an extra 
piece proportional to the fermion hypercharge due to the U(l)-U(l)x mixing. Hence we 
need to re-examine the electroweak precision data using the full couplings as well as taking 
into account the effects due to virtual Z' exchanges. 

As mentioned above, there are two interesting limits with very different phenomenologies 
depending on whether Mx is greater or less than Myy, which correspond to high or low scale 
new physics. 

1. M x > M w 

From Eq. (fTTi) we can see that f|<e and the mixing can be ignored. The most stringent 
constraints come from electroweak precision measurements at the Z-pole. As for M0ller 
scattering, the exchange of the Z' gives negligible effect because its large mass and very 
small coupling to the electrons (see Eq. (125]) ). The dominant effect comes from the modified 
Zee coupling, and we find 



where the Thomson limit value of the weak angle, = 0.23867, is used. Unless otherwise 
specified, we will use this value of for all numerical presentations below. We point out 
here that this result is independent of M' . Also if s v and s e have the same sign, 5lr is 
negative. 

2. M x <C M w 




2 ZC-qS^Sw S e 
W 



stv + 1 



-(4 + 15.80s^ + 31.85c,s„s e ) , 



(26) 



In this case we have a light Z' whose mass is essentially Mx (see Eq. ( TT2|) ). and the 
mixing angle rj can be large compare to e. Since the Z' is light, 5lr is dominated by the 
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FIG. 5: Constraint from the light Z' contribution to Alr- The contours are that for the SLAC E158 
lower limit (solid), 5lr = —0.337, and the projected JLAB lower limit (dashed), 5lr = —0.0136. 
All limits are given at 95% CL. The allowed region is in between the respective contours. No 
parameter space exists for 5Alr > 0. 



photon- Z' interference term: 



1 - As 2 Ml 



(4 - 15.80c 2 v s 2 e + 31.85c„s„s 6 ) 



M 2 Z 

Mir 



(27) 



where u^ ,R = — [s v Q z ' (e) + ^c v sws € Y e y From the SLAC E158 measurement of Am, we 
have at 95% CL 

f M \ 2 

s 2 q - lh.8Q C y e + 31.85c v s v s e G [-4.06, 1.47] x 1(T 5 f J . (28) 

If we use the projected JLAB limits from above, the constraint interval on the RHS is change 
to [—0.16,0.16] in the same unit. Now from Eqs. (1111) and (112H . s e can be written in terms 
of M z > and as 

(29) 



c 2 W M 2 , 
M 2 W 



Using this, we plot in Figs[5]the allowed (s v , M z >) region as given by Eq. fl28|) . Note that for 
the RHS of Eq. (128]) greater than zero, there is no real solution, and so no allowed parameter 
space. 
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The mixing parameters r/ and e are also constrained by low energy observables such as 
the anomalous magnetic moment of the muon a M . The new physics contribution from a light 



extra Z' to a M is found to be 



111 



5a, 



G^ml Ml 



[s 2 {As% - 2s 2 w - 1) - 2(3s 2 w - 2)s w s rj c ri s e + s 2 w c 2 v s 2 ] 



6x/2tt 2 Mf, 
-1.94 x 1(T 9 (4 - 0.19c 2 s 2 - 1.00c v s v s e ) 

-1.94 x 10~V 1 1 



Mf 

mi: 



M 2 - r 2 M 2 



/If 4 s 2 



M^(s w - 0.19) + 0.19MJ, 



Mf 
M 2 ' 



(30) 



where we have used Eq. ( 129]) in the last line. We see that for Mz' -C Mw, the curly bracket 
is positive definite. Thus a 5a^ arising from a light Z' is always negative. This means that 
if confirmed, the current BNL result of a 3.4<r deviation from SM expectation 24J: 



cxp 



(2.95 ±0.88) x lQ- 



Cai) 



cannot be due to the new physics arising from the light Z' alone. 



E. The MCRS Model 



Recently, the warped five- dimensional (5D) RS model endowed with a bulk gauge group 



25j - which we cal 



26 



27 



the MCRS model - has been 



28|. The gauge hierarchy problem 



SU(3) C x SU(2) L x SU(2) R x U(l) x 
actively studied framework for flavor physics 
is solved when the warp factor is taken to be kirr c pa 37 16|. This sets the scale for the 
new physics to vtikk ~ 4 TeV where the lowest Kaluza- Klein (KK) excitation first arises 4 . 
The bulk SU(2)r factor protects the p parameter in the effective 4D theory from excessive 
corrections coming from KK excitations, and is broken by orbifold boundary conditions on 
the UV boundary (brane); this is a 5D generalization of the case studied in Sec. IIII Bl A 
brief summary of the model is given in Appendix |AJ which sets our notations. 

In the MCRS model, 5lr arise due to the interactions of KK excitations of the neutral 
gauge bosons and the charged leptons. Contributions come either from direct virtual ex- 
changes of neutral gauge KK modes, or through the Zee coupling been modified due to 
mixings of SM modes with KK modes. The latter involve both fermions and gauge bosons, 



Note this is at the higher end of the LHC discovery limits. 
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and is similar to the shadow Z' case discussed above. More details of the vertex correction 
due to KK mixings can be found in Appendix [A] 

A generic feature of the extra dimensional models is that gauge-// couplings are not 
flavor diagonal in the mass eigenbasis. To calculate the contributions to 5lr in the mass 
eigenbasis requires the knowledge of both the LH and RH rotation matrices, which are 
determined by the mass matrix of the charged leptons and their localization in the extra 



dimension (see Ref. 26j). In Table [TTJ we list the contributions pertinent to the case of 
polarized M0ller scattering 5 . We give results for five typical charged lepton configurations 



Config. 


Direct KK exchan 
LL (1(T 5 ) RR (1(T 6 ) 6 


g£(io- 4 ) 


LL (10- 4 ) 


Zee correction 
RR (10- 5 ) 


$£f (10- 3 ) 


Total 

Wio- 3 ) 


I 


1.98574 


1.51504 


8.09458 


-1.07684 


-7.98263 


-2.45877 


-1.64931 


II 


1.98956 


1.51484 


8.11154 


-1.07827 


-7.98167 


-2.47224 


-1.66109 


III 


1.99197 


1.51050 


8.12409 


-1.07920 


-7.96023 


-2.49935 


-1.68694 


IV 


1.99188 


1.51065 


8.12365 


-1.07916 


-7.96172 


-2.49767 


-1.68530 


V 


1.99049 


1.51483 


8.11566 


-1.07858 


-7.98162 


-2.47504 


-1.66347 



TABLE II: New physics contributions to the polarized M0ller scattering in the MCRS model for 
five typical charged lepton configurations that give rise to experimentally observed charged lepton 
masses, and are admissible under current LFV constraints. 

that lead to charged lepton mass matrices compatible with the current lepton flavor violation 
(LFV) constraints, and give rise to experimentally observed charged lepton masses. Details 
of the charged lepton configurations and the particular realization of the admissible charged 
lepton mass matrices used to calculate the new physics contributions in Table [Til are given 
in Appendix [B] Note that the direct KK LL (RR) contribution corresponds to the Wilson 
coefficient c' LL (crr) scaled by a factor (2\^2Gf)~ 1 - 

As can be seen from Table HT1 the new physics contribution to Alr from direct exchanges 
of virtual KK gauge modes are subdominant compare to that from the shift in the SM Zee 

5 We do not displayed the tiny imaginary parts, which is < O(10 -21 ), since the real part of the contribution 
dominates in the tree- level Alr. 
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coupling due to KK mixing effects. This is because the gauge KK modes are heavy 6 and 
couplings of the electron to gauge KK modes are small. Note that the LL-type contribution 
is always greater that of the RR-type. This is a consequence of the fact that the LH charge is 
greater than the RH charge, i.e. |Q§(e)| > \Q§ (e)|. Note also the new physics contributions 
are quite stable across the configurations, indicating their robustness with respect to the 
variations in the charged lepton mass matrix. The MCRS values of 5lr at vhkk ~ 4 TeV 
are well within the projected JLAB limits of \5lr\ < 0.0136. 



F. Doubly Charged Scalars 

Doubly charged scalars P ±± are motivated by Type II seesaw neutrino mass generation 
mechanisms. They can be either SU(2)l singlets or triplets, and they carry two units 
of lepton number. The contribution of P ±=t to M0ller scattering depends on only two 
parameters, viz. its mass, Mp, and its coupling to electrons, y ee . For example, if there is 
only coupling to RH electrons, the effective contact interaction is 

l ^(^e R )(e R e R ) = ^(e^Re^Re) , (32) 

leading to 

|2 



>LR 



\y. 



• (33) 

2^2(1-4406^! 

For models where the coupling is to LH electrons, 8l R is similarly given with the RR Wilson 
coefficient changed to that of the LL type. The important difference is that Slr is negative 
here. We plot the Slr arising from doubly charged scalars coupling only to RH electrons 
in Fig. El as well as the relevant constraints from the SLAC E158 and the projected JLAB 
1 ppb measurement. 



IV. CONCLUSIONS 



We have examined in detail the implications of a very high precision polarized M0ller 
scattering experiment on new physics beyond the SM, which is well known to be sensitive to 
extra Z' and doubly charged scalars in the few TeV mass range. We have shown how such 

6 As required by electroweak precision measurements and very stringent quark sector flavor constraints, 
m K K > 4 TeV 



17 



$LR 
1.000 - 
0.500 - 

o.ioo ; 

0.050 - 



E158 



0-010 . ' JLAB 
0.005 - " " 



0.001 



0.15 0.20 0.30 0.50 0.70 1.00 



\Ve 



FIG. 6: The ratio, 5lr, as a function of \y ee \- The solid, dashed, and dotted lines denote cases where 
A = 1, 2, 4 TeV respectively. The horizontal lines marks the SLAC E158 upper limit, 5lr = 0.122, 
and the projected JLAB upper limit, 5lr = 0.0136. All limits are given at 95% CL. 

a measurement can be particularly useful if such a state is found at the LHC. If a doubly 
charged scalar were discovered at the LHC, a measurement of the 5lr will determine its 
Yukawa couplings to electrons, and the sign of Slr will determine whether it couples to LH 
or RH electrons (negative for the former, positive for the latter). Such scalars also carry 
two units of lepton number, hence measurements of 5lr will in turn have impact on the 
predicted rates of neutrinoless double beta decays of nuclei jl2|], thus providing invaluable 
information on the origin of neutrino masses. Similarly, if an extra Z' were discovered, 
measuring 5lr will shed light on its origin. We have given detailed calculations for various 
classes of models including the Eq GUT model, left-right symmetric model, and two examples 
of non-anomalous U(l). 

For high scale new physics, we found if the new particle is of the more exotic type such 
as a hidden/shadow Z' or the pure KK Z' in warped extra-dimension RS scenarios, 5lr 
measures the change in the effective Zee couplings as compare to the SM, and it is not 
sensitive to the mass of these new gauge bosons. In the scenario where the new physics is of 
a very low scale (<^C v), and a Z' with a small mass exists that couples very weakly to SM 
fermions, a 5lr measurement gives even more valuable information not available from high 
energy hadron collider experiments. We have shown that the SLAC E158 experiment can 
already tightly constrains the mixings of this type of Z' with the SM Z . The smallness of the 
mixings allowed conforms with string-theoretic expectations, as well as phenomenological 
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determinations 



111 ]. On the other hand, the positive deviation from the SM expectation in 



the measured value of muon g — 2 cannot be accommodated by such a light shadow Z' since 
it gives a negative contribution. To the best of our knowledge there is no cleaner or better 
measurements on couplings of a hidden/shadow Z' to electrons. These are valuable inputs 
to dark matter model building employing a hidden or shadow sector. 

It is clear that an extra Z' or doubly charged scalars with masses in the 1 to 2 TeV 
range can be found at the LHC relatively easily. We have shown that Slr can provide 
cross checks on these measurements as well as independent information on the couplings to 
electrons of different chiralities. If the masses are > 4 TeV, discovery at the LHC will take 
longer, and may have to await an intensity upgrade to the Super LHC (SLHC). Then it is 
even more important to do a high precision 8lr measurement. The current reach of new 
physics of a 10 ppb measurement achieved at SLAC is around 1 to 2 TeV. Given the same 
couplings, a 1% measurement will probe new physics at a scale A ~ 3 — 6 TeV, and so is 
nicely complementary to the SLHC. This is true in particular for the MCRS model, as the 
quark sector flavor bounds already constrained the masses of the lowest KK modes to be in 
at or above this range. 
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APPENDIX A: RS MODEL WITH SU(2) L x SU(2) R x U(l) B -l BULK SYMMETRY 

We describe briefly in this appendix the basic set-up of the MCRS model to establish 
notations relevant for studying flavor physics. A more detailed description can be found in, 



e.g. Ref. |25|. 



The MCRS mode is formulated on a slice of AdS$ space specified by the metric 

ds 2 = G AB dx A dx B = e' 2 ^ rj^dx^dx" - r 2 c d<j) 2 , (Al) 
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where cr((f>) = kr c \(f>\, rj^ = diag(l, —1, —1, —1), k is the AdS$ curvature, and — n < <fi < ir. 
The theory is compactified on an S\/(Z 2 x Z' 2 ) orbifold, with r c the radius of the compactified 
fifth dimension, and the orbifold fixed points at = and <fi = n correspond to the UV 
(Planck) and IR (TeV) branes respectively. To solve the hierarchy problem, one takes 
knr c pa 37. The warped down scale is defined to be k = ke~ knTc . Note that k sets the scale 
of the first KK gauge boson maSS, Tflgauge 

pa 2. A5k, which determines the scale of the new 

KK physics. 

The MCRS model has a bulk gauge group SU(3) C x SU(2) L x SU(2) R x U(l) x under 
which the IR brane-localized Higgs field and transforms as (1, 2, 2) . The SM fermions are 
embedded into SU{2) L x SU{2) R x U{l)x via the five- dimensional (5D) bulk Dirac spinors 

Qi (u, L[+ , +] \ +]\ Di /**[-,+] 

\di L [+, +] / \diR [-, +]) \ d iR [+, +] 

L l= h l+ ' +] ), B , = f^|-' + M, !=1 , 2 ,3, (A2, 
yea [+, +] J \e iR [-, +] J 

where Qi transforms as (2, l)i/6, Ui and A as (1, 2)i/ 6 , Li as (2, l)-i/2, and Ei as (1, 2)_i/ 2 . 

The parity assignment ± denote the boundary conditions applied to the spinors on the 

[UV, IR] brane, with + (— ) being the Neumann (Dirichlet) boundary conditions. Only 

fields with the [+,+] parity contain zero- modes that do not vanish on the brane. These 

survive in the low energy spectrum of the 4D effective theory, and are identified as the SM 

fields. Note that since we are only interested in the charged leptons in this work, we need 

only the RH charged leptons, and it is not necessary to have the doubling in the lepton 

sector as in the quark sector. 

A given 5D bulk fermion field, can be KK expanded as 

e 

/r c 7r 

where subscripts L and R label the chirality, and the KK modes /£ R are normalized accord- 
ing to 

l[d<t>ftMf?M = 5 rnn. (A4) 

The KK-mode profiles are obtained from solving the equations of motion. For zero-modes, 
we have 
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3a/2 00 



VljiM) = ^ (A3) 



n=0 



where the c-parameter is determined by the bulk Dirac mass parameter, m = ck, and the 
upper (lower) sign applies to the LH (RH) label. Depending on the orbifold parity of the 
fermion, one of the chiralities is projected out. 

After spontaneous symmetry breaking, the Yukawa interactions localized on the IR brane 
lead to mass terms for the fermions on the IR brane 



S- 



Yuk 



d x 



kr r ir 



Q(x,ir)\%U(x,ir)+Q(x,ir)\ d 5 D(x,ir) + L(x,ir)\ e 5 E(x,n) + h. c. , (A6) 



where vw = 174 GeV is the VEV acquired by the Higgs field, and \^ d ' e are (complex) 
dimensionless 5D Yukawa coupling matrices. For zero-modes, which are identified as the 
SM fermions, this gives the mass matrices for the SM fermions in the 4D effective theory 



kVpTX 



kr r 7i 



(A7) 



where / labels the fermion species. The mass matrices are diagonalized by a bi-unitary 
transformation 



(U f L yMfU f R 



I m{ ^ 

ml 



0m 



/ 



(A8) 
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where m{ are the masses of the SM quarks and leptons. The mass eigenbasis is then defined 



by ip' = Uty, and the CKM matrix given by V( 



CRM 



Because of KK interactions, the couplings of the Z to fermions are shifted from their SM 
values. These shifts are not universal in general, which leads to flavor off-diagonal couplings 
when the fermions are rotated from the weak eigenbasis to the mass eigenbasis: 



£ FC nc D ^Z, { Qz(f') £ K L ab f aL rf bL + Qz(f') E «S faRl'fm 

a,b a,b 



(A9) 



The mass eigenbasis is defined by /' = Wf, and the flavour off-diagonal couplings are given 
by 

Kaf = E( C/ U)-4' R (^^)^ » ( A1 °) 

where = diag(«i, k 2 , k 3 ) are the coupling shifts due to KK interactions in the weak 
eigenbasis. Note that there would be no flavor violations if k is proportional to the identity 
matrix. 
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The coupling shifts from two sources: mixing between the SM Z and neutral KK 

gauge bosons, and that between the SM fermions and their KK excitations. These mixing 
processes are depicted in Fig. [7J In the gauge mixing diagram, X can be either the SM Z 



f <H> 



< H > 



<H> 




f 



fi 



<H> 




FIG. 7: Effective Zff coupling due to gauge and fermion KK mixings 



or the Z ! that arises from the SU(2)r, while in the fermion mixing diagram, i, j, and a are 
generation indices. Note that the Yukawa interaction mixes LH (RH) SM fermions with the 
RH (LH) KK fermions. One can see from the diagrams in Fig. [7J that depend on the 
overlap of fermion wavef unctions, which are controlled by how the fermions are localized in 
the extra dimension. Details of the calculation and the full expressions of the contributions 



due to gauge and fermion KK mixings can be found in Ref. 



26|]. 



APPENDIX B: TYPICAL ADMISSIBLE CHARGED LEPTON MASS MATRI- 
CES 

Parameterizing the complex 5D Yukawa couplings as As^- = p^e^, admissible mass 
matrices of the forms given by Eq. (IAT[) are found with pij and 0^ randomly generated in the 
intervals (0, 2) and [— n, tt) respectively. In the following, for each of the five charged lepton 
configurations determined by the c-parameter sets cl and ce, we list for a typical solution 
that passes all current LFV constraints the mass eigenvalues, the magnitude of the complex 
mass matrix, \M\, of the solution, and the magnitude and phase of the corresponding 5D 
Yukawa coupling matrix 7 . All numerical values are given to four significant figures (except 
for the mass eigenvalues, which required more accuracy to distinguish m T ). The mass 
eigenvalues agree with the charged masses at 1 TeV found in Ref. [29| to within one standard 

7 We need not list arg M because it is the same as <fi, as can be seen from Eq. (|A7[) . 
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deviations quoted. 
• Configuration I: 



c L = {0.6265, 0.5562, 0.5002} , c E = {-0.6911, -0.5806, -0.5192} 
{m e , m M , m T } = {4.95902 • HT 4 , 0.104688, 1.77976} 



( 0.0003248 0.004970 0.04274^ 



|M| 



V 



( 0.6598 0.2733 0.4445^ 



P = 



0.3272 0.2883 0.4624 
0.6022 0.4513 0.5351 



0.001416 0.04607 0.3908 

0.009891 0.2738 1.717 

/ 

,0 = 



0.2714 -0.8414 -3.055 
-1.509 2.910 -1.787 
-1.294 1.547 



Configuration II: 



c L = {0.6474, 0.5802, 0.5172} , c E = {-0.6714, -0.5566, -0.5026} 
{m e , m M , m T } = {4.95902 • 10" 4 , 0.104688, 1.78001} 



IMI = 



A 



/ 



Configuration III: 



7999 


0.4591 


0.5226 


5945 


0.2770 


0.6256 


4253 


0.3419 


0.5346 


III: 







0.0003856 0.008502 0.03564 
0.002473 0.04427 0.3681 
0.009663 0.2984 1.718 

\ Z-2.562 -1.623 

-1.491 -0.6484 -0.05579 
-1.274 2.565 



3.094 



c L = {0.6903, 0.5694, 0.5165} , c E = {-0.6269, -0.5683, -0.5029} 
{m e , m M , m T } = {4.95902 • 10" 4 , 0.104688, 1.77966} 

^0.0006886 0.002653 0.009256^ 

IMI = 



0.0006886 0.002653 0.009256 
0.01672 0.1594 0.5872 
0.03637 0.1366 1.670 
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i 1.375 0.8461 0.5746^ 

0.6662 1.014 0.7273 

yO.3606 0.2163 0.5147 

Configuration IV: 



/ -2.730 -2.449 0.7525^ 







-0.9022 -0.2020 -3.037 
1.814 2.844 



c L = {0.6908, 0.5786, 0.5142} , c E = {-0.6277, -0.5583, -0.5057} 
{m e , m M , m T } = {4.95902 • 10~ 4 , 0.104688, 1.77990} 



( 0.0007459 0.002838 0.01386^ 



IMI 



V 



I 1.556 0.6898 0.9268 

1.009 1.108 0.9381 

^0.3642 0.2838 0.5176 

• Configuration V: 



0.01879 0.1771 0.5452 
0.03762 0.2517 1.669 

\ 



I 1.701 1.106 -1.161 
-1.361 -1.094 2.471 
-2.407 0.6682 



c L = {0.6477, 0.5689, 0.5002} , c E = {-0.6714, -0.5681, -0.5190} 
{m e , m M , m T } = {4.95902 • 10~ 4 , 0.104688, 1.78002} 



IMI 



A 



1 1.268 0.4201 0.3330 
0.1955 0.2316 0.4316 
yO.5200 0.5493 0.5265 j 



0.0006072 0.005543 0.01604 
0.001143 0.03729 0.2536 
0.01666 0.4850 1.696 

\ 

0.04108 -1.505 -1.032 
-1.560 2.667 



/ 

( 0.5492 -2.288 1.186 ^ 
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